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ABSTRACT 

We use the XMM-Newton EPIC-pn instrument in timing mode to extend spectral 
time-lag studies of hard state black hole X-ray binaries into the soft X-ray band. We 
show that variations of the disc blackbody emission substantially lead variations in 
the power-law emission, by tenths of a second on variability time-scales of seconds or 
longer. The large lags cannot be explained by Compton scattering but are consistent 
with time-delays due to viscous propagation of mass accretion fluctuations in the disc. 
However, on time-scales less than a second the disc lags the power-law variations by 
a few ms, consistent with the disc variations being dominated by X-ray heating by 
the power-law, with the short lag corresponding to the light-travel time between the 
power-law emitting region and the disc. Our results indicate that instabilities in the 
accretion disc are responsible for continuum variability on time-scales of seconds or 
longer and probably also on shorter time-scales. 
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1 INTRODUCTION 

Rapid X-ray variability on time-scales of seconds to mil- 
liseconds is a key characteristic of X-ray binaries (XRBs) 
which has been studied almost since their discovery. It was 
suggested quite early on that the rapid variability observed 
in XRBs and on much longer time-scales in active galax- 
ies, was produced by instabilities in the 'standard' accretio n 
disc (|Lightman fc Eardley|[l974l ; IShakura fe Sunvaevlll97fJ) . 
Such a model is attractive because the disc can produce 
variations over a broad range of time-scales, corresponding 
to the range of unst able disc radii, to match well with the 
observed variability (Lvubarskii 1997). 

However, for many years the variability 
has been associated with the still-mysterious 
hard X-ray power-law emitt ing component (e.g. 
Galeev. Rosner. fc Vaianal Il979l: iPoutanen fc Fabian! 



and hard spectral states, the disc-dominated soft states 
show ver y little rapid variab ility (fractional rms ~ 1 per 
cent, e.g. [Belloni et al. I l2005l) . In contrast, rms of several 
tens of per cent is seen in the hard spectral states which 
are energetically dominated by the power-law emission 
thought to be associated with a hot optically-thin flow, 
magnetised corona, or even the base of the persistent 



19991: IChurazov. Gilfanov. fc Revnivtsevl l200lt 

Done. Gierlihski. fc Kubotal 2007 and references therein). 
This is because, although strong variations in the lumi- 
nosity of the disc blackbody component are seen on long 
time-scales corresponding to transitions between the soft 



radio -e mitting jets seen i n the se states (e.g. Zdziarski et al.l 
19981 : iDroulans etTafl l20ld : iMarkoff. Nowak. fc Wilms] 
2005). In addition, the strong QPOs observed during the 
state-transitions of black hole XRBs (BHXRBs), when 
both disc and power-law emission are strong, show hard 
rms-spectra consistent with an origin predominantly in th e 
power-law emitting component |Sobolewska fc Zvckil feOQG). 

In the hard states our understanding of the origin of 
the variability has been limited by a lack of spectral cov- 
erage of the disc, which emits at kT < 0.5 keV and so is 
not covered by the bandpass of detectors such as the Pro- 
portional Counter Array (PCA) on the Rossi X-ray Timing 
Explorer, which is the prime workhorse of X-ray timing stud- 
ies to date. Recently we addressed this problem by studying 
the X-ray variability 'covariance spectraQ of the hard state 
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The covariance spectrum measures the spectral shape of varia- 
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of th e black hole candidate GX 339-4 (jWilkinson fc Uttlevl 
I2009T ) using the soft X-ray response of the EPIC-pn CCD 
detector on board the XMM-Newton satellite. In that work, 
we found that the disc blackbody emission does vary with 
an rms amplitude of tens of per cent and that the black- 
body variability amplitude is greater relative to the corre- 
lated power-law variations on longer time-scales (> Is). 
We used these results to argue that, although on time-scales 
< 1 s disc blackbody variability could be explained by X- 
ray heating of the disc by the varying power-law emission, on 
longer time-scales the disc is intrinsically variable. A natural 
physical interpretation is that instabilities in the standard 
disc are responsible for driving the observed power-law vari- 
ability on time-scales at least down to a few seconds, perhaps 
through mass-accretion variations which propagate through 
the disc |Lvubarskii|[l997l ; lArevalo fc Uttlevl 120061 ), before 
reaching the power-law emitting region. 

An even stronger argument can be made for disc-driven 
variability if we can establish the causal relationship between 
correlated variations in the disc and power-law components, 
through the detection of X-ray time-lags between variations 
in different energy bands. Previous work to measure time- 
lags in BHXRBs has used proportional counter instruments 
to measure lags between variations at harder energies, abov e 
2 keV (e.g. lMivamoto fc Kitamotolll989l ; rNowak et all 19991 ). 
however we can also do this analysis using the same XMM- 
Newton EPIC-pn timing mode data that we use to study 
the covariance spectrum. In this Letter we use these data 
to present the first time-lag study of the hard state which 
is extended to soft X-rays to cover the disc component. We 
confirm that there is indeed a clear signature in the lags at 
soft energies which shows that variations on time-scales of 
seconds or longer in the disc blackbody emission lead cor- 
related variations in the power-law emission. We also show 
that the sign of the lag changes on time-scales < 1 s, so that 
the disc component lags behind the power-law variations by 
a few milliseconds. These lags are consistent with light-travel 
time-lags expected by disc thermal 'reverberation', caused 
by X-ray heating of the disc by the power-law which dom- 
inates over the intrinsic disc blackbody variability on these 
short time-scales. 



2 OBSERVATIONS AND DATA REDUCTION 



We analysed EPIC-pn (jStriider et alJ l200ll ) timing mode 
data from the 2004 March 16-19 XMM-Newton observa- 
tions of GX-339 in a stable hard state. The data were re- 
duced in the standard manner using SAS 10.0.0, processing 
the raw data products using the SAS tools epsplitproc 
and epfast, taking account of background flaring and ex- 
tracting only events with RAWX from columns 31 to 45 
using the SAS tool evselect. Although the tool epfast 
was used on the events file, the Charge Transfer Inefficiency 
(CTI) corrections were insufficient to remove features at 1.8 
and 2.2 keV (Si and Au edges) in the 2004 spectrum. Al- 
though the source is relatively bright, previous analyses have 
shown that the EPIC-pn data are not significantly piled up 
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Figure 1. Top panel: GX 339—4 2004 hard state lag versus fre- 
quency dependence for the medium-soft and hard-medium energy 
bands. Positive values correspond to harder photons lagging softer 
photons. The inset shows the lags on a linear y-axis in the fre- 
quency range where the medium-soft lags drop below the hard- 
medium lags (the x-scale is matched to the scale on the main 
figure). Bottom panel: Power spectra for the soft, medium and 
hard energy bands. 



ijDone fc Diaz Trigdlioiol ) . The SAS tools arfgen and rm- 
fgen were used to generate the ancillary response file (ARF) 
and redistribution matrix file (RMF) for the 15 columns 
used for source extraction. Background was not accounted 
for as it is very small compared to the source and in any 
case background variations are uncorrelated with those of 
the chosen reference energy bands and so do not contribute 
to the lags or covariance spectra. The total exposure time 
for the combined 2004 data was 155 ks. 



3 ANALYSIS AND RESULTS 

3.1 Power and time-lags versus frequency 

We used the event lists, together with Good Time Inter- 
vals (GTIs) excluding telemetry dropouts, to extract light 
curves in three geometrically-spaced energy bands (soft 0.5- 
0.9 keV, medium 2-3 keV and hard 6-9 keVfl in order to 
measure the cross-spectra and determine lags between bands 



tions which are correlated with a given reference band, and so is 
ideal to search for correlat ed disc and power-law variability, see 
IWilkinson fc Uttlevl d2009h for details. 



2 There is a factor ~ 3 separation in average event energy be- 
tween each band, with values of 750 eV, 2420 eV and 7140 eV for 
soft, medium and hard bands respectively. 



© 2010 RAS, MNRAS 000, [UU 



Connection between disc and power-law variability 3 



which are separated by equal steps in log(_E). Cross-spectra 
for pairs of these energy band s were measured in the stan- 
dard way ijNowak et alJll99Sh for contiguous segments of 
data of length 59 s (the segment-length necessitated exclud- 
ing data containing short gaps due to telemetry drop-outs). 
The cross-spectra were averaged over segments and in geo- 
metrically spaced frequency bins with frequency separation 
v — > 1.15^. We measured phase lags from the argument of 
the average cross-spectrum in each bin and divided by 27r 
times the average frequency of the bin to give the time-lag. 
The resulting lag-frequency spectra are shown in Fig. [1] to- 
gether with the power spectra of the measured energy bands. 

From the lag-frequency spectra, we note two striking 
differences between the lag behaviour in the soft band and 
the behaviour at harder energies which has already been 
well-studied by RXTE. First, although we chose the energy 
bands to be separated by roughly-equal steps in log(J5), at 
low frequencies the medium-soft lags are a factor ~ 6 times 
larger than the hard-medium lags. Therefore, by extending 
the measurement of lags into the soft band we have un- 
covered a significant increase in the lags compared with 
the values expected from extrapolation of the log-linear 
scaling with energ y observed by RXTE above 3 keV (e.g. 
iNowak et al.l [19991 ) . which would give similar lags between 
pairs of bands separated by the same \og(E). It is also inter- 
esting to note that although at low frequencies the medium- 
soft and hard-medium lag-frequency relations show simi- 
lar power-law shapes (time-lag r oc v~ ' 7 , as previously 
observed), at frequencies approaching 1 Hz the medium- 
soft lags start to drop rapidly, to the point (above 3 Hz) 
where they drop below the hard-medium lags. This drop 
in medium-soft lags seems to correspond roughly to a drop 
in the soft band power relative to oth er bands at frequen- 
cies ap proaching and above ~ 1 Hz. In I Wilkinson fc Uttlevl 
(2009) we used the covariance spectra measured on differ- 
ent time-scales to associate this drop in soft band power at 
high frequencies with a drop in the amplitude of variability 
intrinsic to the blackbody-emitting disc, so we might expect 
the disc to play a similar role in explaining the medium-soft 
lag behaviour. 

3.2 Lag versus energy and covariance spectra 

To better understand the origin of the lag behaviour associ- 
ated with the soft band, we carried out a Fourier-frequency- 
resolved study of lag versus energy. To optimise signal-to- 
noise we measured the cross-spectra for many narrow en- 
ergy bins with respect to a fixed broad reference band cov- 
ering the range 0.54-10.08 keV (the reference energies are 
selected to match the start and end of energy bins). In order 
to remove the correlated Poisson noise part of the variabil- 
ity from the cross-spectrum we must remove the energy bin 
light curve from the reference band light curve before mea- 
suring the cross-spectrum. Strictly-speaking this means that 
the reference band is slightly different for each energy bin, 
but due to the narrow bin widths, the effects on the rela- 
tive lags are negligible compared to the erro rs (see also the 
discussion by IZoghbi, Uttlev. fe Fabiar]|201Ch . 

After measuring the cross-spectrum for each channel 
with respect to the reference band, we can determine the 
relative lag versus energy for a selected frequency-range 
by averaging the cross-spectrum over that range (e.g. see 



INowak et alii 19991 ; I Kotov. Churazov. fc Gilfanovll200lh . We 
can use the frequency-averaged cross-spectra to determine 
phase lags and hence time-lags in the usual way. We can also 
use the averaged cross-spectrum to make frequency-resolved 
covariance spectra, which can be obtained by dividing the 
square-root of the modulus-squared of the cross-spectrum by 
the rms of the reference band, in a manner analogous to the 
way the time-domain covarian ce spectrum is measured us ing 
the cross-correlation function (I Wilkinson fc UttlevH2009l ). 

Fig. [5] shows the frequency-resolved lag-energy spectra 
for four frequency ranges, each corresponding to an approxi- 
mately factor 4 increase in frequency starting with 0.034 Hz 
and ending at 8 Hz. Covariance spectra for the same fre- 
quency ranges are shown as insets and are plotted as ratio 
plots with respect to a simple absorbed power-law model 
(see Figure caption for more details). 

The frequency-resolved lag-energy spectra show clearly 
that the soft excess emission associated with the disc black- 
body is responsible for the complex soft lag behaviour seen 
in Fig. [1] The increased amplitude of the low-frequency 
medium-soft lags results from a significant downturn in the 
lag-energy relation, away from the log-linear law which ap- 
plies above 2 keV. The covariance spectra show that this 
downturn is associated with energies where variable disc 
blackbody emission starts to become significant. Therefore, 
the disc blackbody variations are significantly leading the 
variations of the power-law emission which dominate at 
harder energies. 

The covariance spectra show that at high frequencies 
the disc blackbody variability is still present but drops in 
amplitude. The sa me result was found in th e time-domain 
covariance spectra (I Wilkinson fc Uttlevll2009h , our interpre- 
tation being that at these frequencies, intrinsic disc vari- 
ability is small compared to blackbody variations driven by 
X-ray heating of the disc by the power-law emission. This 
interpretation appears to be confirmed by the lag-energy 
spectra. The downturn in low-frequency lag-energy spectra 
becomes an upturn at high frequencies, so that now the soft 
photons lag the medium-energy photons (which explains the 
sharp change in medium-soft lags seen above 1 Hz in Fig.[T]). 
This behaviour can be simply explained if the mechanism 
for disc blackbody variability is switching from being dom- 
inated by variable viscous heating intrinsic to the disc on 
time-scales below 1 s, to being dominated by X-ray heating 
by the illuminating power-law on shorter time-scales. The 
high-frequency soft lags, which correspond to tens of Rg 
light-crossing time, can be associated with the light travel 
time from the power-law emitting regions to the blackbody- 
emitting regions of the disc. 



3.3 Lags in other hard state systems 

To show that other hard state systems show the same low- 
frequency lag behaviour as GX 339—4, we plot in Fig. |3]lag 
versus energy for the 0.125-0.5 Hz range for EPIC-pn timing 
mode data from XMM-Newton observations of hard states 
in GX 339-4, Swift J1753. 5-0127 and Cyg X-l as well as 
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Figure 2. Lag-energy spectra for four temporal frequency ranges. The lags are all measured relative to the same 0.54-10.08 keV reference 
band and plotted so that energies with more-positive values of lag are lagging behind energies with less-positive values of lag. The insets 
show the covariance spectra for the same energy range and frequencies, plotted as a ratio to a power-law with photon index fixed at 
r = 1.55, absorbed by a 6 X 10 21 cm -2 neutral Galactic column and fitted to the 3—10 keV range. Note that the scales for the covariance 
plots are the same, showing that the switch in soft lag behaviour at high frequencies corresponds to a drop in amplitude of disc variability 
relative to the power-law variability. 



for GX 339—4 in a lower-flux state than observed in 
The lags at higher frequencies cannot be compared with the 
GX 339—4 2004 data due to larger error bars caused by the 
shorter exposure times of these observations (40 ks, 16 ks 
and 31 ks for Swift J1753. 5-0127, Cyg X-l and GX 339-4 
2009 respectively). However, in the 0.125-0.5 Hz range the 
other systems, as well as the GX 339—4 2009 data, show a 
downturn in lag versus energy in soft X-rays which is similar 
to that seen in the GX 339-4 2004 data. 



4 DISCUSSION 

By extending the measurement of energy-dependent lags to 
soft energies, we have shown that on time-scales of seconds 
or longer, disc blackbody variations lead variations in the 
power- law component. The natural implication of this result 
is that the disc variations drive the power-law variations, so 
we now consider the physical origin of this causal connec- 
tion. The medium-soft lags of 0.1 s or longer seen in vari- 
ations on time-scales of seconds correspond to light-travel 
times of thousands of gravitational radii. Coronal Comp- 
tonisation models have been invoked to expl ain the log- 
linear form of the lag s seen at harder energies (|Pavnelll980l ; 
iKazanas fc Huall 19991 ) , but they struggle to explain even the 
more moderate hard-medium lags seen at low frequencies 
without invoking an unfeasibly large Comptonising region 

3 Swift J1753. 5-0127 was observed on 2006 March 24, GX 339-4 
and Cyg X-l were observed on 2009 March 26 and 2009 December 
2 respectively. 
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Figure 3. Lag-energy spectra of 0.125—0.5 Hz variations for sev- 
eral hard state systems. 

(|Nowak et all I1999T I. A recent more energetically-feasible 
model invokes Comptonisation b y the persistent hard state 
jet to produce the har d lags (|Reig. Kvlafis. fc Gianniosl 
120031 ; iKvlafis et afl feoOS'). In this model, seed photons are 
fed into the base of the jet (e.g. from the disc), and the 
anisotropy of Compton scattering along the jet axis means 
that upscattered photons can traverse large distances in the 
jet before escaping to the observer, to produce large lags rel- 
ative to the photons undergoing fewer upscatterings. How- 
ever, in such a model we would expect the medium-soft 
lags to be reduced relative to the hard-medium lags, not 
increased as observed here, since the direct seed photons 
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from the disc at soft energies will dilute the log-linear lags 
introduced by the upscattered continuum. The medium-soft 
lag would only increase relative to the hard-medium lag if 
the medium-soft lag represents the actual travel time from 
the seed photon source to the scattering electron population, 
i.e. the seed photons must first interact with the jet not at 
its base but thousands of gravitational radii above the disc. 
This geometry raises significant problems for explaining the 
energetics and variability of the emission (since the jet at 
that height would subtend only a small solid angle as seen 
from the disc), not to mention the shape and strength of the 
iron Ka line, which due to geometric and beaming effects 
would be much narrower and weaker than observed, if the 
disc were illuminated by th e jet from such a great height 
(e.g. iMarkoff fe Nowakll2004h . 

A likely explanation of the large medium-soft lags 
observed on time-scales of seconds is that the lags 
are associated with the generation of and propaga- 
tion of accretion fluctuations within the in ner regions 
of the disc before they reach the corona dLvubarskiil 
1997 : IKotov. Churazov; fc Gilfanovll200ll : lArevalo fc Uttlevl 



20061 '). In fact, the observed lag behaviour on these and sub- 



second time-scales is exactly that expected from our ear- 
lier disc- variability interpre tation of the covariance spectrum 
I Wilkinson fc Uttlevl [20091 ). In the present work, we have 
used causal information to greatly strengthen that inter- 
pretation, so that accretion instabilities can be firmly iden- 
tified as the source of X-ray variability, at least on time- 
scales of seconds. Furthermore, the unstable accretion flow 
is now shown to be the standard disc, w ith all the attendant 
phys ical implications (see discussion in IWilkinson fc Uttlevl 
120091 ). and not a hot optically thin flow. 

The disc propagation model can explain the medium- 
soft lags, however the hard-medium lags - produced where 
the disc blackbody does not contribute to the X-ray spec- 
trum - still require explanation. These lags may still be 
related to Compton upscattering of seed photons as the 
accretion fluctuations reach the base of the jet, e.g. in 
a hybrid model of disc propa gation and the model of 
Rci E . Kvlafis. fc Gianniosl (|2003h . However, given the com- 
mon i/ -0 ' 7 frequency-dependence of the low- frequency lags 
which seems to be independent of the energy-bands chosen, 
it seems likely that the lags at hard energies are directly re- 
lated to the same underlying propagation mechanism which 
produces the medium-soft lags. For example, lags at harder 
energies can be simply produced if the mass-accretion fluc- 
tuations in the disc preferentially generate softer and then 
harder power-law emission as they propagate inwards. Such 
effects could be produced if there is a soft coronal component 
abov e the disc while the hard emission is produced centrally 
(e.g. IKotov. Churazov. fc Gilfanovl l200ll ; lArevalo fc Uttlevl 
l2006h . 

The switch in lag behaviour on short time-scales sug- 
gests that the disc lags the power-law variations by a few 
ms on these time-scales, consistent with the light-travel lags 
expected from a power-law component separated from the 
disc by only tens of gravitational radii at most. This lag sig- 
nature represents the first evidence in BHXRBs for a disc 
thermal reverberation lag. This is the time-lag due to the 
light travel-time from the central power-law emitting region 
to the disc where the hard X-rays are reprocessed into black- 
body emission. In order to use these lags to map the disc, 



e.g. to constrain the inner radius, we would need to make 
assumptions about the geometry of the power-law emitting 
region and the emissivity profile of the disc blackbody. This 
modelling effort is beyond the scope of this work. 

Finally, it is important to note that although the lags 
observed for variations on time-scales of less than a second 
appear to be caused by X-ray heating of the disc, this does 
not necessarily imply that the variability on these shorter 
time-scales is not also generated in the disc. An accretion 
instability at small disc rad ii will modulate only di sc emis- 
sion inside that radius (e.g. lArevalo fc Uttlevll2006h . so the 
intrinsic variability in disc emission will be small and likely 
to be dominated by X-ray heating effects when the mass 
fluctuations reach the central power-law emitting region. 
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